The Brownian gyrator: a minimal heat engine on the nano-scale 
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A Brownian particle moving in the vicinity of a generic potential minimum under the influence 
of dissipation and thermal noise from two different heat baths is shown to act as a minimal heat 
engine, generating a systematic torque onto the physical object at the origin of the potential and 
an opposite torque onto the medium generating the dissipation. 
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Introduction and summary: The theory of heat engines 
is one of the main roots of modern thermodynamics and 
statistical physics. Recently, there has been a consid- 
erable renewed interest in the long standing problem of 
determining the fundamental efficiency limit of a heat en- 
gine at maximum power [l|. Another exciting new per- 
spective is the conceptual design of a cooling device on 
the nano- or even single-molecule scale by inverting^a 
heat engine which is powered by Brownian motion [3]. 
The original and still paradigmatic setup of a heat en- 
gine consists of two heat baths in contact with a cycli- 
cally working "engine" , generating work in the form of 
a torque. Here, we put forward the smallest and most 
primitive such engine one may think of: a single particle, 
gyrating around a generic potential energy minimum un- 
der the influence of friction and thermal noise forces from 
two simultaneously acting heat baths. Three typical ex- 
amples are indicated in Fig. 1. While the particle itself 
is way too small to store or disburse any notable amount 
of angular momentum, it acts as a kind of catalyst. The 
particle in fact generates, with the help of the disequilib- 
rium between the two baths, a systematic average torque 
onto the physical object at the origin of the potential and 
an opposite torque of the same magnitude by way of the 
dissipation mechanism onto one or both heat baths. 

Such Brownian gyrators are minimal heat engines in 
so far as they are acting essentially like their macro- 
scopic counterparts, but the "engine" itself needs not 
be anything more than a structureless particle perform- 
ing Brownian motion. Due to this simplicity they are 
in principle readily realizable by nano-technological or 
even single-molecule techniques. Potential applications 
are manifold and obvious: wrapping and unwrapping of 
DNA and other polymers fsl , driving wheels or screws of 
nano-devices Q and synthetic molecular motors [^l^stir- 
ring and mixing in micro- and nano-fluidic devices [5| , to 
name but a few. Though each single engine is "weak" , 
when acting in parallel, a case naturally arising in the 
context of colloidal particles [1, 0| or magnetic fluxons 
Q , the result may even be a macroscopic torque 0] . 

Model: Focusing on the simplest case, we consider the 
motion of a point particle with two spatial degrees of 
freedom x :— {xi^x^) in a static potential J7(x). We 



assume that the potential has a minimum at the origin 
X = and that large excursions are sufficiently rare to 
admit a parabolic approximation of the form 

f/(x)=^|y^ y,:=^0,,(aK, (1) 

i=l fc=l 

where 0{a) is a 2 x 2 orthogonal matrix with elements 
Oil = O22 — cos a, O12 = —O21 = sin a, describing a 
rotation in the plane by an angle a. The transformed 
coordinates yi thus correspond to the "principal axes" 
of the parabolic potential in ([T]) and ui, M2 > are the 
corresponding "principal curvatures" . Clearly, Eq. ([T]) 
represents the generic form of a potential minimum in 2 
dimensions. Furthermore, one has ui ^ U2 in the generic 
case, i.e. unless the system is rotationally symmetric 
about the origin. 

To complete our heat engine, we need two heat baths 
which act onto the particle without resulting in a to- 
tal equilibrium system. In the most common case, the 
heat engine is alternatingly brought into contact with 
two baths at two different temperatures. Since this re- 
quires a quite complicated machinery in practice, here we 
rather focus on the case that the particle is permanently 
in contact with both baths. In the following, we first dis- 
cuss in detail the conceptually simplest theoretical model 
and only afterwards turn to the experimental realizations 
and more general system classes. 

In the simplest case, at least one of the two heat baths 
interacts with the particle along a preferential direction, 
which can be identified with the xi-axis without loss of 
generality. Generically, this direction is not related in any 
particular way to the "principal axes" of the potential in 
, and hence their relative angle a is not bound to take 
any particular value a priori. 

Whereas one heat bath thus solely influences the par- 
ticle motion along the xi-axis, the second one may either 
be of isotropic character or acting only on X2- For the 
sake of simplicity, we first focus on the latter case, see 
Fig . l(a,b). Modeling the thermal bath effects as usual 
[l3) [m by Gaussian white noise and a concomitant dissi- 
pation proportional to the velocity, and neglecting inertia 
effects, we arrive at the following overdamped Langevin 
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FIG. 1: Schematic sketch of various types of Brownian gy- 
rators acting as heat engines. The black dot represents the 
Brownian particle, moving in the a:i-a:2-plane. The contour 
lines indicate a typical parabolic potential ([1} with principal 
a:xes y\ and j/2. (a)-(c) illustrate different realizations of the 
two heat baths, (a): Both heat baths act on the charged par- 
ticle by way of black body radiation at different temperatures 
T\ and T2, irradiating along the x\- and a;2-axes, respectively. 
The dissipation mechanism is provided by radiation damping 
into the vacuum, (b) : "Electrical heat baths" , realized by two 
resistors at different temperatures T\ and Tb, coupled to the 
charged particle by means of two plate condensers. Each of 
them transfers the random voltage fluctuations of one resis- 
tor to the particle along a preferential direction and gives rise 
to dissipation via the resistor when the particle moves and 
hence induces a current in the electrical circuit. Replacing 
the condenser plates by Helmholtz coils gives rise to "mag- 
netic baths" interacting with e.g. a paramagnetic particle [9|. 
Replacing the condenser plates by piezo elements gives rise 
to "acusto- mechanical" baths [ij. (c): Only one heat bath 
(with temperature T') is of the anisotropic type as in (a) and 
(b). The second heat bath (with temperature T) consists of 
the usual fluid environment of the Brownian particle. 



equations for the particle dynamics in the plane: 

ac/(x(i)) 



+ V2^^J^^^rit) , 1= 1,2 . (2) 



Here, fcs is Boltzmann's constant, T,; is the tempera- 
ture of the i-ih bath, ^i{t) are independent, (5-correlated 
Gaussian noises, and the coupling strength between par- 
ticle and bath i is quantified by the friction coefficient r]i 

MM- 

Solution: We first discuss in some more detail the 
forces and torques connected with the dynamics 
Denoting by the unit vector along the i-th coor- 
dinate axis, the three relevant forces are: the dissi- 
pation frjit) ~J2i=i^iViii{t): the potential force 
iu{t) ■— — V?7(x(i)), and the fluctuation force fj(i) 
J2i=i ^iV^^Vi^BTi £.i{t) ■ Hence, ^ is tantamount to the 
force balance f^(i) -I- iu{t) + i^it) = 0. On the aver- 
age over many realizations of the noise we thus obtain 
{{^{t)) = and hence {iu{t)) — — (fi)(i))i corresponding 
to the following elementary physics: Since the particle 
momentum is by definition considered as negligible in 
the overdamped limit, the force exerted by the potential 
U is compensated on the average by the friction forces. 
A similar consideration applies to the three torques of 
the form f x x: The thermal fluctuations do not give rise 
to any systematic torque, {i^{t) x x(i)) = 0, and hence 
the torque of modulus M{t) and direction 63 := ei x 62 
which the particle exerts on the potential U (or the physi- 
cal object at the origin of that potential) is on the average 
exactly equal to the opposite torque —M{t)e3 which the 
particle exerts via the friction forces on the thermal envi- 
ronment (or the physical objects containing the baths): 

(fuit) X x(i)) = -(f^(t) X x(t)) = M{t)e3 . (3) 

Next we turn to the Fokkcr-Planck-equation [l^, [ll| 
equivalent to 



dP{x,t) 
dt 
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where /'(x, t) is the probability density to find the par- 
ticle at position x at time t and J = (Ji, J2) the corre- 
sponding probability current density with components 



1 dU{x) keTi d 
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The Fokker-Planck-equation (I4I5P is complemented by 
natural boundary conditions P(x, i) —^ and Ji{yi,t) — > 
for Xi ±00. Given P(x, t), the torque modulus M{t) 
from ([3]) readily follows according to 



M{i) = J P(x,i) 



xi- 



dU{x.) at/(x) 



dx2 



X2- 



dxi 



dxi dx2 (6) 



After initial transients have died out, the system ap- 
proaches a unique, steady probability density as t — > 00 
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[Tol. [Til. [T3| . For the parabolic potential ([T]), this unique 
steady state solution of the Fokker-Planck-equation (|4l5p 
can be obtained in closed analytical form. Since the cal- 
culations are straightforward but rather tedious and the 
expressions quite bulky and not very illuminating, they 
are not explicitly given here. Rather we immediately 
present the resulting torque ^ in the steady state, read- 
ing 

_ ksiTi - T2) {ui - U2) s\n2a 
U1+U2- {ui - U2) cos 2a 

Note that for symmetry reasons, any value of the po- 
tential force fc/(0 occurs with the same probability as 
its inverse in the steady state. On the average, we thus 
have {fu{t)) = 0, implying as usual that the resulting 
torque ([7]) remains unchanged for any other choice of the 
reference rotation axis in ([3]). 

Discussion: The closed, general expression ([7]) for the 
average torque of a Brownian gyrator in the steady state 
is the main result of our paper. For Ti = T2 we are 
dealing with an equilibrium system in ^ and hence 
the average torque must vanish due to the second law 
of thermodynamics. If ui ~ U2 the potential H]) is 
rotationally symmetric and hence there cannot be any 
preferential direction of rotation and the torque must 
vanish. If sin2Q! = then the principal axes of the 
parabolic potential agree with the directions along which 
the two heat baths in ^ are acting, hence a net torque 
is again ruled out by symmetry. In any other case a fi- 
nite torque ([7]) results (the denominator is always positive 
since u^, 77^ > 0). A particularly simple behavior arises 
for equal coupling strengths 771 = ri2- In this case, the 
maximal torque is reached at a = 7r/2 + wk. In general, 
the optimal angles will be slightly different. Typically, 
the maximal torque is roughly given by the difference be- 
tween the thermal energies associated with the two baths, 
fcs Ti — fcs r2 , indicating that the Brownian gyrator trans- 
forms its random motion into a systematic torque quite 
effectively. We, however, remark that speaking about ef- 
ficiencies in the usual sense is not possible as long as one 
does not know the resulting rotation speed of the ther- 
mal baths relatively to the "carrier" of the potential J7, 
which is is beyond the scope of the general model ^ . 

The basic physical origin of the preferential rotation 
of the Brownian particle in one direction can be most 
easily understood in the limit that one temperature van- 
ishes, say T2 = in ([2]). Now, let us assume the particle 
has reached (in whatever way) a position on the xi-axis, 
i.e. X2{t) = 0. In the generic case that the xi-axis does 
not coincide with a principal axis of the parabolic po- 
tential in ([1]), there will be a non- vanishing deterministic 
force —dU{xi{t),{))/dx2 proportional to xi{t) acting on 
the particle along the a;2-direction. Since the noise along 
this direction vanishes in we can conclude that the 
particle is able to cross the xi-axis only in one direction 
for all positive a;i(t)-values and in the opposite direction 
for all negative xi(i)-values. In other words, the particle 
rotates around the origin in a preferential direction. It is 



plausible that qualitatively an analogous behavior is ex- 
pected also for finite T2 (different form Ti), though the 
details will be more complicated. 

Some of the above basic physical principles governing 
the behavior of a Brownian gyrator are similar to those 
governing so-called thermal ratchets and Brownian mo- 
tors [ll|. Yet, to the best of our knowledge, neither of 
those ratchet systems are immediately comparable to the 
setup treated here, nor are the general concepts in the 
context of ratchet effects of any help to gain easier or 
deeper insight in the present case. 

Experimental realizations: The setup from Fig. 1(a) 
involving anisotropic black body radiation is of consid- 
erable conceptual interest. The main practical problems 
are the weak coupling of a charged particle to the electro- 
magnetic irradiation and the vacuum, and that the model 
(|2|) itself is a very crude description of the real system. 
Yet, the basic concept may well be of relevance for var- 
ious dynamical processes in astrophysics, space physics, 
and laboratory plasmas, see pisj and references therein. 
More easy to realize in the lab is the setup from Fig. 1(b), 
whose heat baths consist of simple resistors at different 
temperatures. If the so generated thermal fluctuations 
are still too weak, an electronic amplification is straight- 
forward [3]. The shortcoming of this setup is that a 
conversion of the torque into a relative rotation between 
potential and baths is not desirable, since that would 
change the angle a in ([1]). 

Most attractive from the experimental viewpoint seems 
to be the setup from Fig. 1(c): One heat bath is given 
- as usual in the context of Brownian motion - by a 
fluid environment of the particle without any kind of 
anisotropy. For later use, we denote its temperature by T 
and the Stokes friction coefficient by 77. They both appear 
as usual [13, [HI in both components of the 2-dimensional 
dynamics ([2]). Only the second bath continues to emit its 
thermal fiuctuations along a preferred direction onto the 
particle, say along the xi-dynamics. Paradigmatic exam- 
ples are [1^ the (almost) black body irradiation from the 
sun or some analogous artificial device in the lab or any 
kind of anisotropic (almost) white noise in the original 
mechano-acoustical sense, emitted e.g. by a loudspeaker 
in the experimental work [l4| . Denoting its temperature 
and dissipation coefficient by T' and rj' , respectively, the 
resulting 2-dimensional dynamics is again of the over- 
damped Langevin type and can be readily brought into 
the form ([2]) by means of the identifications 772 — rj, 
7^2 = r, 771 = 77 + T2 = {vT + iT')/{rj + 77'). The 
resulting torque ([7]) thus takes the form 

ksiT' - T) (ui - U2) sin2a 
M = ^ . (8) 

U1+U2- 2^p;7 (mi - U2) COS 2a 

Typically, the coupling strength rj to the fiuid will be very 
much larger than the dissipation coefficient 77' due to the 
anisotropic second heat bath. In order that the relevant 
strength 

9 keT'Tj' (9) 
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of the concomitant anisotropic fluctuations is non- 
negligible in spite of the small 77'-value, the temperature 
T' must be very much larger than T . Then the resulting 
torque from ([S]) simplifies in very good approximation to 

M = ^^^i^sin2a. (10) 

Tj Ul + U2 

Here, any non-zero value of M indicates that the particle 
is transferring torque from the potential U to the dis- 
sipative mechanism, while its own angular momentum 
always remains negligible. The effect of that pair of op- 
posite torques will be to generate rotations of the phys- 
ical object carrying the potential and the fluid around 
the particle into opposite directions, just in the way any 
"engine" is commonly supposed to operate. 

The experimental realization of the potential U is pos- 
sible in many straightforward ways, e.g. by means of 
electro- or magnetostatic forces dielectrophoretic ef- 
fects (including light forces as exploited in optical traps) 
0, pinning centers of fluxons Q, etc. 

Outlook: We close with some generalizations and per- 
spectives. First of all, basically the same behavior is 
expected when working in 3 rather than 2 dimensions. 

More challenging is to better understand the role of the 
anisotropy of at least one of the heat baths with the main 
goal of possibly abandon this condition. The first pur- 
pose of this anisotropy is of mainly technical character. 
Namely, within the usual modeling of thermal fluctua- 
tions as Gaussian white noise and the concomitant dis- 
sipation proportional to the instantaneous velocity, see 
two isotropic baths are effectively equivalent to one 



single equilibrium environment with properly adapted ef- 
fective friction and temperature. Hence the anisotropy is 
indispensable in order to obtain a non-equilibrium sys- 
tem within this standard modeling of the thermal baths. 
However, there are in principle many possibilities - some 
of mainly conceptual interest, others of practical rele- 
vance but mathematically more difficult to handle - to 
model a thermal bath in a different way, e.g. by means 
of correlated Gaussian noise and concomitant memory 
friction Jllj. Two of these baths at different tempera- 
tures can no longer be mathematically transformed to 
one single equilibrium bath and hence in this regard the 
anisotropy is no longer needed. A second more funda- 
mental role of the anisotropy is to break the symmetry 
between gyrating clock- and counterclock-wise. Clearly, 
breaking this symmetry is an indispensable pre-requisite 
of making the Brownian gyrator work. Hence, in the 
presence of two isotropic heat baths, this symmetry must 
be broken in some other way. The most straightforward 
possibility is via the potential U, e.g. by keeping terms 
up to cubic order in the expansion ([T|) of the potential 
about its minimum. 

An interesting extension of the present work will be to 
explore the collective phenomena due to many interact- 
ing Brownian gyrators, in particular the similarities and 
differences as compared to collective effects of rotating 
molecular motors in membranes [161] and rotating mag- 
netic discs confined to a two-dimensional interface [iTj . 
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